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ABSTRACT 



An equation is derived to calculate the dynamics of relativistic magnetized plasma which 
decelerates by sweeping up matter from the ISM. Reduction to the non-radiative and radiative 
^ . regimes is demonstrated for the general case of a collimated plasmoid whose area increases as a 

r| I power of distance x from the explosion, and in the specific case of a blast-wave geometry where 

the area increases quadratically with x. An equation for the evolution of the electron momentum 
distribution function in the comoving fluid frame is used to calculate the observed synchrotron 
radiation spectrum. The central uncertainty involves the mechanism for transfering energy 
from the nonthermal protons to the electrons. The simplest prescription is to assume that a 
l/^ , fixed fraction of the comoving proton power is instantaneously transformed into a power-law, 

C^ ' "shock" -like electron distribution function. This permits an analytic solution for the case of a 

relativistic plasmoid with a constant internal magnetic field. Effects of parameter changes are 
f~^ , presented. Breaks in the temporal behavior of the primary burst emission and long wavelength 

0^ ' afterglows occur on two time scales for external ISM density distributions next oc x~^ . The first, 

as emphasized by Meszaros & Rees, represents the time when the outflow sweeps up « A4th/ro 
Mh| of material from the ISM, where A^th is the mass in the outflow ejecta and Fq is the initial 

Q ' Lorentz factor of the plasmoid. The second represents the time when synchrotron cooling begins 

strongly to regulate the number of nonthemal electrons that are producing radiation which is 
observed at a given energy. 

Results are applied to GRB and blazar variability. The slopes of the time profiles of the 
X-ray and optical light cuves of the afterglows of GRB 970228 and GRB 970508 are explained by 
1^^ ' injection of hard electron spectra with indices s < 2 in a deceleration regime near the radiative 

j^ ■ limit. We also propose that blazar flares are due to the transformation of the directed kinetic 

energy of the plasmoid when interacting with the external medium, as would occur if relativistic 
outflows pass through clouds orbiting the central supermassive black hole. 



Subject headings: galaxies: active — galaxies: jets — gamma-ray bursts — gamma rays: theory 
— radiation mechanisms: nonthermal 



1. Introduction 

CGRO BATSE and Beppo-SAX observations of GRBs have renewed interest in relativistic blast-wave 
models. The BATSE observations (Meegan et al. 1992; Fishman & Meegan 1995) indicate that GRBs 
originate from cosmological distances and therefore involve impulsive releases of energy exceeding 10^^ ergs 
(assuming isotropic emission). The Beppo-SAX observations of X-ray afterglows from GRB 970228 (Costa 
et al. 1997) and GRB 970508 (Piro et al. 1997) led to the first optical counterpart identifications of GRB 



sources (van Paradijs et al. 1997; Djorgovski et al. 1997), culminating in the detection of absorption lines 
from the optical counterpart to GRB 970508, showing that it lies at redshift z > 0.835 (Metzger et al. 
1997). The X-ray emission from GRB 970228 decays as a power-law ex t~^^^ '^'^'^ between w 50 s and one 
week following the GRB (Costa et al. 1997), where tobs is the observer time. The optical counterpart of 
GRB 970228 decays ex iotV^^"^""® between 0.6 days and 189 days following the burst, though it exhibits 
apparent fluctuations from a simple power-law (Galama et al. 1997; Garcia et al. 1997). 

The X-ray afterglow of GRB 970508 displays a peak at about 2 days following the burst, with an 
underlying power law decline oc t^^^^ from the end of the main phase of the burst at about 30 sec to « 6.1 
days following the GRB (Piro et al. 1997). The optical flux of GRB 970508 appears to be roughly constant 
between « 6 hrs and 1 day, whereupon it subsequently rises and likewise peaks at « 2 days after the burst 
(Djorgovski et al. 1997; Sahu et al. 1997 and references therein). The optical emission thereafter declines 
ex t;7bs^*^°°'' (Garcia et al. 1997). The 4.86 GHz radio emission from GRB 970508 was first detected « 6 
days after the GRB, is rapidly variable (evidently due to interstellar scintillation; Goodman 1997) though 
with an overall flat time profile until « 60 days following the burst, after which it begins to monotonically 
decline (Frail et al. 1997). 

Because of these observations, much interest has been focused on the hydrodynamical study of 
relativistic blast waves by Blandford & McKee (1976), which has been applied to GRBs by Rees, Meszaros, 
& coworkers in a series of papers (e.g., Meszaros & Rees 1993, 1997; Meszaros, Rees, & Papathanassiou 
1994; Wijers, Rees, & Meszaros 1997; Panaitescu & Meszaros 1997) devoted to understanding the prompt 
GRB emision and the long wavelength afterglows. Vietri (1997a, b) used the results of Blandford & McKee 
in the radiative regime to explain the delayed X-rays, and Waxman (1997a, b) applied the blast wave model 
to the long wavelength afterglows. Meszaros, Rees, & Wijers (1997) have derived power-law time profiles 
for delayed emission in different regimes assuming analytic models for the ISM density distribution. 

Here we reframe the hydrodynamical approach by deriving an equation for plasmoid dynamics in 
terms of the comoving particle distribution functions. This permits the radiation physics to be treated in 
the comoving frame, and the plasmoid dynamics to be calculated self-consistently. This approach can be 
applied, with appropriate modifications, to temporal studies of blazar emissions which are also thought to 
result from relativistic plasma outflows. Blazar flares, which vary on time scales as short as hours to days 
in the 1 MeV - 1 TeV regime (e.g., Michelson et al. 1994; Catanese et al. 1997), are in this picture due to 
the reconversion of directed kinetic energy of the outflowing plasmoid by interactions with the ISM (see 
Hartman et al. 1997 and Shrader & Wehrle 1997 for recent reviews of blazar observations). We suggest that 
the passage of the relativistic plasma outflow through a cloud-like density enhancement could initiate such 
a flare. 

In §2, we derive an equation for plasmoid dynamics in terms of the evolving comoving particle 
distribution functions. Reduction to well-known forms in the non-radiative and radiative regimes is 
performed in §3. Prescriptions for the energy transfer from nonthermal protons to electrons, which get us 
past the difficult particle acceleration and plasma physics which can be implemented in later studies, are 
proposed in §4. In §5, we treat the optically thin synchrotron radiation process, and effects of parameter 
changes are presented in §6. The results are examined analytically in §7, and application to GRBs and 
blazars is made in §8. We summarize in §9. A detailed numerical treatment is presented in a companion 
paper (Chiang & Dermer 1997b). 



2. Plasmoid Dynamics 

We consider the system shown in Fig. 1. A plasmoid, consisting of a volume of magnetized plasma, 
is located on the x-axis at location x and moves with Lorentz factor r(a:;). The comoving volume of the 
plasmoid is V{x) = A(x)L(x), where A{x) is the area of the plasmoid at x and L(x) is the mean width of 
the plasmoid as measured in the comoving frame. The plasmoid travels through an external medium with 
density next (a;) which is assumed to be fully ionized electron-proton plasma. Consequently the plasmoid 
captures an electron-proton pair of the ISM gas if the proton Larmor radius is much smaller than the 
characteristic size scales {L and A^'^) of the plasmoid. Except for chance collisions, neutral particles would 
pass through the plasmoid unless they are ionized by interactions with the plasmoid. 

The plasmoid is assumed to consist initially of thermal lepton-proton plasma which was ejected from 
the GRB event. In our notation, nip =- 1.5 x lO^'' ergs and rrig = 8.2 x 10^^ ergs, so that the mass of the 
thermal particles in the plasmoid is just A^th — fnpNth{l + ath), where iVth is the total number of protons 
in the plasmoid ejecta. The factor oth = (1 + 2zth)"T^e/'7ip corrects for the mass of electrons and pairs, 
the latter of which are assumed to be present in the ratio zth = A^th,-i-/-^th, where A'th,+ is the number 
of thermal positrons. Here we assume that k^T^/nif^ <C 1 and k^Tp/mp <C 1, so that the thermal kinetic 
energy is a small fraction of the rest mass energy of the particles. 

We furthermore adopt the notation that Npi^\ stands for the proton (electron) momentum distribution 
function Npi^\ [p: x{t)\ differential in dimensionless momentum p = j3^ (not to be confused with the subscript 
p which stands for proton), and that N' ,-. stands for the spatial gradient dNp(^f,-^[p]x{t)\/ dx of the proton 
(electron) momentum distribution function, where t is the time in the comoving frame. If the nonthermal 
particles are isotropically distributed in the comoving frame (see discussion in Dermer & Schlickciser 1993), 
then the relativistic mass A^nt of nonthermal particles is given by 



Mnt ^mp dp J {Np + antNe) , (1) 

Jo 

where ant(p; 2;) = [1 + 2znt(p; x)]me/mp is the additional relativistic mass in the form of nonthermal leptons, 

and Znt{p) is the momentum-dependent pair ratio in the nonthermal Icpton distribution. 

The x-componcnt of the relativistic momentum of the plasmoid at location x, denoted Ilx{x), is given 
by P{Mth + A^nt), where P = BY and Be is the speed of the plasmoid. At x -|- 6x, the total x-component 
of the momentum of the plasmoid, including the swept-up particles and radiated photons, is given by 

n^(x + Sx) = mp{P + P'6x){{Nth + A^th'5a;)(l + oth + a[i^Sx)+ 

f°° 6x f°° 
+ dp-f[Np + N'pSx + (a„t + a'^^Sx) [N^ + N'Jx)] } / dp {%Np + %Ne). (2) 

Jo C Jq 

The last integral in this equation refers to the momentum carried by radiated photons, and is related to the 
internal particle distribution functions through the comoving frame energy-loss rates mpjp and mgje of the 
protons and electrons, respectively. This formalism allows for particle creation and annihilation, although 
in this paper we assume for simplicity that ajj^ = a'^^ = 0. We furthermore set N^.^ = so that no particles 
are added to the thermal pool. This does not impose any loss of generality if the nonthermal particles are 
followed to nonrelativistic energies. 

Expanding equation (2) in powers of 6x, equating to H-xi^) by momentum conservation and retaining 
terms to first order in dx, we obtain the equation 

P' _ /q°° dp j{N; + a„t^^) - (cP)-^ l^ dp {%Np + %N,) 



P Nthil + ath) + j;^ dp -f{Np + antNe) 



(3) 



The proton and electron injection function due to the plasmoid sweeping up particles from the ISM is 
given by 

7V;,^ = n,,t^J(p-P)=.iV;,^. (4) 

Equation (4) says that protons and electrons are injected into the comoving plasma with momenta p = BY 
at a rate determined by the ISM density next (2;) and area A{x) of the plasmoid. When radiation is 
important, the functions N' and N'^ in equations (2) and (3) contain terms in addition to the sweep-up 
functions (4) due to the evolution of the comoving particle distribution functions. By multiplying the 
particle continuity equation by 7 and integrating over ^7, using the relation dx = BTcdt, we find that 

d7 7{JV(7)+ Q^ } = J^ d7[7iV(7)-77V(7)] = Brcy^ ^7 7^sw , (5) 

noting that ['J'jN{'y)]\^ — 0. Consequently the numerator in equation (3) contains only the contribution 
from the sweep- up function (4), and we arrive at the result 



P' ^ /o°° dp 7(JV;,^ + antiV^sw) 

P iVth(l + ath)+/o°°rfp7(^^P + ant^e) 



(6) 
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Note that we have neglected the energy content in magnetic field in equation (3) by assuming that it is 
small in comparison with that of the thermal and nonthermal particle energies. The importance of the 
magnetic energy term can be checked once the plasmoid volume is specified, and can be important if 
equipartition is assumed between the magnetic field and nonthermal particle energy densities. For the 
applications presented here, however, we assume that this term is negligible. 



3. Asymptotic Behavior of r(x) 

The relation between the comoving time t and spatial coordinate x is given by 

The terms xq and ig in equation (7) stand for the initial location and comoving time when the parameter 
values are specified. The nonthermal power from swept-up protons and leptons is probably dominated by 
the proton component in the ISM environment surrounding a GRB source, though a substantial pair content 
could be found in the vicinity of supermassive black holes which power blazars. Assuming for simplicity 
that the swept-up matter is pair free, we can therefore write that the power injected in nothermal particles 
per unit comoving time is given hy E = nip L dp 7 Ap.sw = "^p^rc/„ dp 7 N' ^^. Using equation (4), we 
therefore sec that the kinetic energy injected per unit time is 

i^kc = mpBiV^ - r)cnextA (8) 

(Blandford & McKee 1976), neglecting the power in swept-up electrons. Vietri (1997a, b) has used this 
equation to predict GRB afterglow behavior assuming that the bulk of the nonthermal particle energy is 
radiated on a time scale short compared to the time scale for the evolution of F, i.e., that the plasmoid is 
in the radiative regime. 

When particle radiation is negligible, 7p = 7e = 0, and the evolution of the particle distribution 
function is entire due to the sweep-up nonthermal particles. Under these conditions, equation (6) reduces to 

_P^ ^ no^t{x)A{x)r{x) 

P " iVth + Jo di next(x)A(x)F(£) ■ 



Equation (9) describes the dynamics of the plasnioid in the non-radiative regime u. If the particle kinetic 
energy is promptly radiated away, leaving behind only the rest mass of the initial and swept-up particles - 
which is dominated by the mass of the protons - then equation (6) reduces to 

P Ntb + /g^ dx nc^t{x)A{x) 
This equation describes the dynamics of the plasnioid in the radiative regime. 



3.1. Non- Radiative Regime 

Asymptotic forms for plasnioid dynamics in the non-radiative and radiative regime are simply derived 
for a relativistic (F ^ 1) plasmoid. We parameterize the plasmoid area by the expression 

A = A^{-y^Anxlh{-y , (11) 

Xq Xo 

where /;, = (551/47r is a coUimation (or beaming) factor and 5Vl is the solid angle into which the plasmoid is 
directed. Equation (11) reduces to an uncoUimated blast wave when j = 2 and fb = ^- We parameterize 
the spatial dependence of the external density distribution by the expression 

next(xO=no(-)-", (12) 

Xo 

and assume that the plasmoid Lorentz factor follows the power-law behavior 

^(:r) = ^o(-)-^ (13) 

Xo 

where Fq is the initial plasnioid Lorentz factor. 

In the early free-expansion stage, the plasmoid is ballistic so that F(a;) = Fq. Transition to power-law 
behavior in both the non-radiative and radiative regimes occurs when the relativistic mass Mbvj^o of the 
swept-up particles equals the original mass mpNth in the burst ejecta. This occurs when 

dx Hc^tix) A{x)T{x) = Fo / dx ncxt{x)A{x) > iVth = 7; , (14) 

Jo f o"^p 

where the equality on the right-hand-side relates the initial energy Eo ~ lO^'^E^o ergs of the GRB explosion 
or blazar flare to A^th- One finds that 

.0 = [^^^^]^/^ = 3.9 X 10-[l^±i-^] V3 em , (15) 

where F300 = Fq/SOO. This reduces to the result of Rees & Meszaros (1992) when j = 2, 77 = 0, and /(, = 1. 
Substituting equations (11)-(13) into equation (9), one obtains 

_^_^i^J_l^l±l^ (16) 

IX X 



^This limit is also commonly referred to as the adiabatic regime. We avoid use of this term which can be confused with the 
early phases of the fireball expansion where adiabatic losses dominate the energy evolution of particles in the comoving frame. 
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so that equation (9) has a power-law solution in the non-radiative regime provided 

9'^9a = ^ • (17) 

Note that the index ga in the non-radiative regime equals 3/2 for the blast-wave case with uniform external 
density (77 = 0). The power-law behavior persists until the plasmoid reaches location Xf defined by 
T(xf) ~ 1. For the asymptotic form (13), we have 

Xf = rj/^xo . (18) 



3.2. Radiative Regime 

The plasmoid Lorcntz factor in the fully radiative regime follows power-law behavior for xq <^ x <^ Xf, 
where Xf is the location at which the mass of the swept-up matter equals the initial rest mass of 
particles in the GRB ejecta. When A^th S> L dx ncxt{S:)A{x), we can simply see from equation (10) that 
-dr/r^ oc x^-'^dx, so that 

9^9r=i + l-V = '29a ■ (19) 

Alternately, we can solve equation (10) directly when the integral in the denominator on the right-hand-side 
of equation (8) can be neglected. For a blast-wave geometry in a uniform external medium, gr = 3 
(Blandford & McKee 1976; see also Rees & Meszaros 1992; Meszaros et al. 1997). The power-law behavior 
in the radiative regime ends at the location Xf determined by the condition 

fXf 
Nth = / dS: noxt(x)A(x) . (20) 

If Xf ^ xq, this implies 

-_l = ( ^09r )V.. . (21) 

a;o TampnoAoXa 

The power-law behavior ends when r(a;/) « 1. Using equations (18) and (19), we therefore find that xq 
is also given by equation (15) in the radiative regime. Figure 2 shows a numerical solution of T(x) along 
with the asymptotic forms in the radiative and non-radiative regimes for our standard parameter values 
Tsoo = 1, J^o = 1 cm^"^, i?5o = 1, j = 2, 77 = 0, and ft — I (see Table 1). 



4. Energy Transfer from Protons to Electrons 

The bulk of the energy injected into a relativistic plasmoid which sweeps up matter from the ISM is 
carried by nonthermal protons, yet in all likelihood the dominant radiation process in GRB afterglows is 
nonthermal electron synchrotron radiation (sec, e.g., Waxman 1997a; Tavani 1997; Katz & Piran 1997; and 
discussion in Chiang & Dermer 1997a). It is also generally thought that the bulk of the radio-through-optical 
blazar emission is nonthermal electron synchrotron radiation. A mechanism transfers the energy from the 
proton to the electron population, for example, through shocks or the excitation of plasma wave turbulence 
by the protons which gets channeled into the electrons through gyrorcsonant coupling. We avoid the 
complicated microphysical details of the energy transfer process by offering two simple prescriptions for the 
transfer of energy from the protons to the electrons. 



In the first prescription, we assume that a fraction S,pe of the kinetic energy power of the swept-up 
nonthermal protons in the comoving frame is instantaneously transformed into a "shock" -like momentum 
power-law electron distribution with index s. Hence 

4 = ipe{x)E = m,K{x) f ""^' dp (7 - l)p-* = m,K{x)f. (22) 

"'Pmin 

Naive shock acceleration theory predicts s > 2 and s — 2 for strong shocks in nonrelativistic gases (e.g. 
Jones & Ellison 1991). Consideration of relativistic gases and inclusion of nonlinear effects can produce 
harder power law injection with s > 1 (e.g. Ellison, Jones, & Reynolds 1990), and stochastic gyroresonant 
acceleration can produce extremely hard power-law and even quasi-monoenergetic spectra (see Schlickeiser, 
Campeanu, & Lerche 1993). In equation (22), the normalization coefficient K = ^pi,E / frue, where 

./ - / rfp (7 - 1)P-' + r^^ dp (7 - l)p-' = ^^7^^ + ^T ~ ^ ' (23) 

provided Pmin < 1 and pmax > 1- If Pmin > 1, then 



2(3 -s) 



»2-^ - p27S 
J" c^ J^max ^min ('24~l 

A second prescription for the energy transfer from protons to electrons is 

E, = ^^, (25) 

where 

/•OO 

^p.kc = mp / dp (7 - l)A^p . (26) 

Jo 

Note that the energy transfer time scale r^e and the total proton kinetic energy ii^p.kc are x- (or equivalently 

t-) dependent quantities. The second prescription is probably a better representation of the actual physics 

of the situation but affords less analytic development. The limit of complete, prompt energy transfer from 

the protons to the electrons corresponds to S,pe -^ 1 and Tpe ^ in the first and second prescriptions, 

respectively. 

The electron injection function at time ti in the comoving frame for the first prescription is therefore 

Ne{Pi,ti)^ K p^''Q{pi]p^ir,,Praa.yi) = -E^—E. 2^^ p- ^ 9 (p . ; p^. ^ ^ p^^^ ) ^ (27) 

where Pi is the injection momentum and 6(p;Pmin,Pmax) is the Heaviside function defined by Q{x; xi, X2) = 1 
for xi < a; < a:;2 and Q{x; xi, X2) = otherwise. 

The evolution of the electron energy through synchrotron radiation is governed by the well-known 
equation 

4 H^ 
-7 = -ccttuhP^ = vp^ , UH = z 9 ' (28) 

3 STTTOeC 

assuming an isotropic electron pitch-angle distribution. We also suppose that the magnetic field is randomly 
oriented with mean field strength H in our spectral calculations. In general, H = H{x). The evolution 
of H with X involves subtleties. The assignment of an equipartition strength in GRB studies depends on 
knowledge of the local density which can be obtained through the shock jump conditions (e.g., Waxman 



1997a; Meszaros et al. 1997). The magnetic field could be generated through turbulent field growth, or 
evolve from the original field in the plasmoid (Meszaros et al. 1994), for example, by flux freezing (see 
Chiang & Dermer 1997a). Our major simplification in the present study is to examine the case H — Hq 
and v = Vi^ = otcH^ / ((Sirrai^c^) , with _ffo ^nd i'q independent of x. We are also mostly interested in the 
radiation from rclativistic electrons (p 3> 1), in which case the electron equation of motion — p = v^p^ has 
the solution 

V^v{t)^[vi^+Mt-h)r^ (29) 

(see Dermer & Schlickeiser 1993 for a related treatment; Sturncr, & Schlickeiscr 1997 for a treatment 
of Coulomb, bremsstrahlung, and Compton losses in addition to synchrotron losses; and Gould 1975, 
Goldshmidt & Rephacli 1994, Dermer & Skibo 1997, and Chiang & Dermer 1997a for a consideration of 
adiabatic energy losses of the nonthermal electrons) . Here we treat only the effects of synchrotron losses on 
the evolution of the electron spectra. 

The electron distribution at time t resulting from the superposition of nonthermal electron power-law 
injection functions described by equation (27) over the time interval t^ < ti < t is therefore 

N,{p;t) = j dU Ne{p,U)^p-^ I dU K [p-i - Mt - U)Y-^ , (30) 

Jto J max [to, t- 1^0 ^ (P^^ ^PmL)] 

where we assume that Pmin and Pmax are independent of t. We note that a value for Pmax can be obtained 
in shock acceleration theory by balancing radiative loss and energy gain rates (e.g., Reynolds 1996). 

We now follow the evolution of the nonthermal electron distribution that occurs over the period when 
the plasmoid Lorentz factor is described by the power-law asymptote (13). Integrating equation (7) from 
the onset of the power-law behavior at a;o to location Xi, we find that 

^ = Kt,-to) + l]^/(^+^', (31) 

xo 

where u = {g -\- l)roc/a;o. In the range xq < Xi < Xf, T, A, and rioxt sue described by equations (13), (11), 
and (12), respectively. Thus 

mef Xo 

when r > 1. 

In this regime, equation (27) becomes 

^pe mpcTln oAo 
mefp' 
where 

Isir)^ dy (1 + u;r - Luyrip-' - i^oyy-^ ■ (34) 

Jo 

Here t = t — to and u — [j — 2g — rf)/{\ + g) (the integration variable y ~ t — ti). The simplest case is 

given by s = 2, which corresponds to the strong shock index in the quasilinear regime of shock acceleration 

theory. When s — 2, equation (34) becomes 

( [lo{1 + u)]-^[{l + lotY" ^ {I + ojT - ojyY"], if cq 7^ 0; 
h{T) = { (35) 

V^-'Htt^^)^ if ^0 = 0, 

where y — min[T, v^ (p^^ — Pmax)] and cq = 1 + u = {gr — g)/{g + 1), recalling definition (19). 



JVe(p;r)^ ^^-"r?X'^" /.(r) (33) 



5. Synchrotron Radiation 

The synchrotron emissivity in the comoving frame is treated in the (5-function approximation of Dermer 
& SchHckeiser (1993; eq.[5.16]), which we write as 

e Jo CH 2e^ e en 

The dimensionless cyclotron energy en = H/Hcr, where Her — 4.414 x 10^^ Gauss. For a power-law electron 
distribution, approximation (36) is accurate everywhere except near the endpoints of the synchrotron 
spectrum. Self-absorption of the synchrotron radiation and synchrotron self-Compton emission (see Dermer, 
Sturner, & SchHckeiser 1997 for a treatment of the latter process) are not considered here. 

The relationship between the observer's time element 6tohs and the comoving time element 6t is 
Stohs = St{l + z)/V, where the Doppler factor V = [T{1 — B^ohs)]~^ and arccos/Zobs is the angle between 
the direction of travel of an emission element and the observer's line-of-sight. We consider the case where a 
portion of the plasmoid is directed along the line-of-sight to the observer, so that the received radiation is 
produced primarily by emitting plasma which is directed within a line-of-sight angle arccos/iobs ^ 1/r (see 
Chiang & Dermer 1997a when this is not the case). Thus we can let Stohs ~^ St{l + z)/T. a 

The observed photon energy (in units of me) is given by fobs = ^^/(l + z) —^ re/(l -I- z), where the last 
relation again holds for the emitting plasmoid directed within an angle arccos/iobs <^ 1/r of the observer's 
line-of-sight. The number of photons produced per unit observer's time per unit observed photon energy 
Cobs is therefore given by 

-^syn(eobs,iobs) = -^syn,com[e(eobs), ^(iobs)] , (37) 

noting that the Jacobian is unity. The quantity Ncomi^it) is the differential photon production rate in the 
comoving frame. 

In equation (37), we replace e by (1 -I- z)eobs/r and t by a function of iobs defined through the expression 

s r dt 
tohs - tobs,0 = Tobs = (1 + z) / — ^ . (38) 

J to i [t) 

In the power-law regime for the plasmoid Lorcntz factor, equations (38), (13) and (31) imply 

UJT = {flTohs + ir - 1 , (39) 

where n = ujTo{1 + 2.g)/[(l + z){l + g)] = (2g + l)cr2/[a;o(l + z)] and ci = (1 + .g)/(l + 2g). 



6. Numerical Results 

We may now calculate the observed synchrotron radiation spectrum using equations (33) and (36) in 
equation (37). The power-law solution for F becomes accurate after observer's time 

T%li2g+1) Mo 



^ At a fixed observer's time, the mean value of fj, for a shell which emits uniformly per unit surface and expands with constant 
Lorentz factor is equal to B and 4_B/(3 -I- B^) for photons emitted isotropically in the comoving frame with energy index a 
equal to and 1, respectively. When /Jobs ~ B, T) = F. 
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following the start of the GRB, which represents only a couple seconds for our standard parameters (see 
eq.[15] and Table 1). A value of pmax — 10 7-^0 ^^ used, which is consistent with a maximum synchrotron 
photon energy <^ 25 MeV in the comoving frame of expected from diffusive shock acceleration (de Jager & 
Baring 1997). The solution is no longer valid when F « 1, and we stop our calculations there. 

Figure 3 illustrates the relationship between F/Fq and x/xq as a function of the quantity rJiobs- As 
expected, a highly radiative blast wave in a uniform external medium (i.e., g — 3) slows down most rapidly 
and travels the shortest distance. When 17 < 3, due either to a non-blast wave geometry with j < 2, a 
thinning external medium with 77 > 0, or inefficient radiation of the swept-up nonthermal proton energy, 
the plasmoid decelerates more slowly. 

The spectra at different observing times and the time profiles at different observed photon energies are 
presented in Figures 4-6. The spectra are plotted at observing times tobs = 1, 10, 10^ sec, etc., from top to 
bottom. The spectra show a spectral hardening from photon energy index a = 1/2 at lower photon energies 
to Of ^ 1 at higher photon energies due to the effects of cooling on the spectra for the index s = 2 used in 
the electron momentum spectrum. Spectra measured at late times display a spectral break at lower photon 
energies because the effects of cooling on the lower energy electrons are only felt at late times. A pile-up 
appears at the spectral break, and this pile-up becomes more pronounced at later times and lower photon 
energies as more electrons are swept into the pile-up region. 

Figure 4 shows the spectra and time profiles using all standard parameters (see Table 1). The standard 
value of 10 Gauss for the magnetic field is arbitrarily assigned in order to provide a spectral break at ss 100 
keV - several MeV, as is typical for many GRBs (e.g.. Band et al. 1993). The parameters in Figure 5 are 
the same as those in Figure 4 except that Hq — 1 Gauss. This has the effect of moving the energy of the 
spectral break to higher energies at comparable observing times. Figure 6 shows a calculation using all 
standard parameters except for Fg = 30, hq = 0.01 cm^'^, and a beaming factor /;, = 0.01. This choice 
delays the onset of the asymptotic power-law regime to late times (compare eq.[40]), and is driven by an 
attempt to explain observations of the peaking of the optical and X-ray afterglow emission of GRB 970508 
at « 2 days following the burst. 

The time profiles of Figures 4-6 can be understood qualitatively as follows: The highest photon 
energies are produced by very energetic electrons. In the comoving frame, only the most recently injected 
high-energy electrons have not been degraded to lower energies by synchrotron cooling. The time profile 
of the decay at MeV - GeV energies is therefore determined by the power in recently injected electrons, 
which reflects the time-dependence of the power in the swept-up matter as the plasmoid decelerates. The 
time profiles in this regime follow a power law behavior ex t^^^^, with slope x ~ 10/7, as explained in the 
next section. The electrons which produce synchrotron radiation at lower observed photon energies take 
much longer to cool. Before cooling is important, these electrons accumulate. The time profiles in this 
regime consequently decay more slowly than in the strongly cooled regime. The slope of the time profile 
reflects both the accumulation of newly injected electrons and Doppler effects on the observer. Because the 
plasmoid is decelerating, synchrotron emission at a fixed observing energy is produced increasingly by the 
fewer higher energy electrons. 

The spectra in Figs. 4-6 are presented in the dimensionless form vLi^/me, and reach values of 
10^2 _ 1054 during the first 1-100 seconds of the GRB. The weakest GRBs observed with BATSE are 
triggered at thresholds of « 0.2-1 photon cm^^ s^^ in the 50-300 keV range. For a mean photon energy 
of 100 keV for GRB sources at a distance of 10^^ cm, roughly corresponding to z = 1, we would therefore 
expect vL^/nie fluxes of « (0.2 - 1) x 47r x 10^^ x 0.1/0.511 ^ 5 x 10^^ - 3 x 10^^. Consequently, the 
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energetics can be explained only if a typical GRB involves at least 10^^ ergs or if /b <^ 0.01 and H, ^ ^ 1- 
10 sec. 

There is a characteristic time fobs when synchrotron cooling starts to influence the total number of 
electrons which produce radiation at a given energy. The value of fobs is longer at lower photon energies, 
and this effect causes the X-ray to 7-ray flux ratio to rapidly increase during the first 100 seconds in Figures 
4 and 5. For this model with constant magnetic field, the light curve of the decaying gamma-ray emission 
asymptotically approaches a power-law with slope characteristic of the cooled regime. 



7. Analytic Interpretation 

Equation (39) can be rewritten in terms of a dimensionless time T = 1 -I- ilrobs through the expression 

l+wT = ri. (41) 

From equations (13) and (31), 

r = FoT^'^^^ , (42) 

where C2 — g/i^g -1-1). We also write 

EH FoCH 

Equation (37) becomes, with equation (33) and (36), 

7Vsy„(f,T) ^ g^ £-3/2t-3-/2/^(t) , (44) 

where 

(1 + z)mef 
In the uncooled regime, r <C I'o^ {p^^ ~ Pmax)- Examination of equation (34) shows that when 

P < Pmax, 

IsCT) -^ r"'=°°l'='i = P-^lT'^o-c, _ ^^ .4g) 

Equation (45) applies for all values of cq ^ 0; when cq = 0, /^"cooicd _ tj^^ci InT. The time-dependence of 
the energy spectrum of the synchrotron radiation in the regime where the emission is produced mainly by 
electrons which have not cooled through synchrotron radiation is given by 

The regime where the bulk of the injected electrons has already cooled is characterized by the condition 

T > t^(7^(p"^ - PmD- When p < Pmax, 

/,(T) ^ 7™°'°^ = , ^ \ Tci.(co-i) ^ (48) 

(s - 1)vq 



and the observed synchrotron flux is given by 

N^°°^°'^{£, T) = ^° £:-(s+2)/2 j-C2{s+2)/2 + ci-(co-l) ^ (49) 
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From the above results, we see that the negative of the slopes of the synchrotron flux time profiles 
produced by a decelerating plasnioid in the asymptotic power law regimes iobs 3> il^^ are given by 

Xuncoolcd = -Co • Ci + C2(s + l)/2 -^ C2(s + l)/2 , (50) 

and 

Xcooicd = -Co • ci + ci + C2(s + 2)/2 ^ ci + C2(s + 2)/2 , (51) 

in the uncooled and cooled regimes, respectively. The expressions on the right-hand-side of equations (50) 
and (51) represent the fully radiative regime where g = gr and, therefore, cq — 0. Table 2 shows the derived 
slopes of the flux time proflles for emission in the fully radiative uncooled and cooled regimes When g < gr, 
which happens when the incoming energy is not fully radiated - and which must be the case in view of the 
fact that there are still some nonthermal electrons which are radiating energy - the slopes are less negative 
than given here, so that the proflles decay more slowly. A numerical calculation (Chiang & Dermer 1997b) 
is required to determine the slopes of the time proflles in this regime. In the extreme non-radiative regime, 
where g = ga = grl"^, the amount of radiated energy is negligible, and the slopes of the time profiles are 
harder by co • ci = g/{'2,g + 1) units. The derived slopes in the fully non-radiative regime are also listed in 
Table 2, as well as the difference in slopes between the cooled and uncooled regimes. These differences are 
the same in the radiative and non-radiative limits. 

The transition from the regime where synchrotron cooling plays a negligible role to the regime where 
it plays a major role in determining the number of nonthermal electrons with momentum p occurs at 
comoving time f « (vqp)^^ , provided p <g; Pmax- This also yields the observer time fobs when synchrotron 
cooling starts to strongly regulate the number of electrons which produce synchrotron emission at some 
observed energy. For times f ^ lo~^ or, equivalently, fobs 3> ^~^, we can use equation (41) to write 
f = Lu^^ {VlfohsY^ and equation (43) to write p in terms of the observed photon energy fobs- After a tedious 
derivation, one obtains the result 

. _ (1 + z)l--^^(.9 + 1)^-34331-2C3(4_1 ^ 106)C3 ^(,3^^^_g)^s ,grig50,(i_2c..)/3 .„^ 

300(2.9 +l)i/o^-a ''' -h^o^ ■ ^ ' 

The term 53 = l/(2ci + ca) = (2.g + l)/(3.g + 2). When g = 3, 






Tobs - 8.5 tri. 91 0.636 ( E^ „ ) ' ^*^<^ • (^3) 



When 5 = 1, 



Tobs - '.U rrl.8^0.6r0.025 ^ T^ „ ' ^^'^ ' y^^' 

-"O ^obs^ 300 ^50gr 

The time at which synchrotron cooling effects begin to play an important role are very weakly dependent 
upon all parameters except Hq and the observed photon energy fobs- The accuracy of equation (53) can 
be checked against Figs. 4 and 5 using gr ~ g ~ 3. In Fig. 6, the Meszaros-Rees timescale fi"^ is longer 
than fobs, so that synchrotron cooling effects are always important by the time the plasmoid Lorentz factor 
enters the power-law asymptote regime described by equation (13). 



8. Discussion 
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8.1. Application to Gamma Ray Bursts 

We use our results to make preliminary comparison with observations of delayed afterglow behaviors 
of GRBs noting, however, that a detailed numerical scheme is required to calculate accurately the time 
profiles of the afterglows. It seems likely that the X-ray and optical GRB afterglow emission is primarily 
observed in the regime where synchrotron cooling plays an important role in the nonthermal electron 
content which produces this radiation (see eqs.[53] and [54]). From Table 2, we see that the slopes of the 
time profiles in this "cooled" regime, whether in the non-radiative limit or the radiative limit, are fairly 
narrowly bracketed. Irrespective of the electron injection index s in the range 1 < s < 3, the slopes x lis 
between the values of 11/14 = 0.78 and 23/14 — 1.64. If the optical and X-ray afterglow indeed occur in 
the regime where synchrotron cooling is important, then s = 2a (see eq.[49]). For GRB 970228, a = 0.63 
(e.g. Waxman 1997a), and for GRB 970508, a = 0.65 (Djorgovski et al. 1997), implying s = 1.3 in both 
cases. Such a spectrum can be produced by nonlinear effects in a relativistic shock (Ellison et al. 1990). 
The radiative regime with g — gr is suggested on energetic grounds (Vietri 1997b). From Table 2, we 
see that a relativistic blast wave decelerating in a uniform medium yields temporal slopes in the range 
17/14 = 1.21 <^ X ^ 37/28 = 1.32 when 1 < s < 1.5. These values are in striking accordance with the 
slopes of the measured afterglow behaviors (see §1), particularly in view of an additional hardening of the 
temporal slopes expected from the fact that the radiative regime can never fully be achieved. 

Against the apparent successes of the relativistic blast wave model in explaining the slopes of the 
afterglow time profiles, several difficulties must be mentioned: (1) The characteristic GRB spectral shape 
cannot be reproduced with the model presented here, which gives a 1/2 power break due to incomplete 
synchrotron cooling rather than a more typical break of about one unit from the hard X-rays to the soft 
gamma rays. This problem can be ameliorated by the inclusion of relativistic shock effects which produce 
time-varying equipartition magnetic fields and shock-heated low-energy cutoffs in the electron momentum 
distribution, and will be dealt with elsewhere (Chiang & Dermcr 1997b). (2) The short time-scale variability 
and the large diversity of GRB time profiles are not accounted for here, and probably cannot result from 
variations in the external medium (Fenimore, Madras, & Nayakshin 1996; Sari & Piran 1995), though 
internal shocks may suffice (see Kobayashi, Piran, & Sari 1997). (3) The tendency of the peak in the lyF^ 
spectrum of GRBs to be in the range 0.1-1 MeV - which is probably not due to selection effects (Harris 
& Share 1998) - remains unexplained by this model. (4) Finally, the peaking of the X-ray and optical 
time profiles two days after GRB 970508 is unexpected. Fig. 6 shows extreme parameters which produce 
a delayed peaking of the X-ray and optical ranges, but implies an accompanying long term gamma-ray 
production phase. The X-ray and optical feature at 2 days after GRB 970508 must therefore originate from 
another effect, for example by relativistic beaming due to a collimated plasmoid (Chiang & Dermer 1997a) 
or by a density enhancement which the expanding shell traverses as it decelerates. 



8.2. Application to Blazars 

The most extreme fiaring behavior of blazars occurs in the gamma-ray regime, with luminosities Lf 
exceeding 10^^L4g,fb ergs s~^ and gamma-ray variability taking place on time scales tf <^ iday days, with 
Lis ~ 1 and iday <, 1- The Eddington limit for a IO^MsMq black hole is « lO'^'^Ms ergs s"^ If the jet 
luminosity represents a significant fraction of the Eddington luminosity and is collimated into 0.1 - 1% of 
the full sky, then 1O*M0 black holes could power such flares. This interpretation is also consistent with a 
minimum variability time scale defined by the Schwarzschild radius of the black hole, which is about 15 
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minutes when Mg ^ 1. 

Interactions of the jet plasma with clouds in the external medium could provide a mechanism for 
transforming the kinetic energy of the jet power into the emitted radiation. Certain conditions have to be 
met for this interpretation to be valid. In particular, the clouds have to be (i) small enough to be in accord 
with the observed variability time scale, (ii) dense and massive enough to afford efficient conversion of the 
plasmoid energy into radiation, (iii) close enough to the central black hole so that a cloud can pass through 
the volume traversed by the jet on a time scale of order days, (iv) large enough to occult a large portion of 
the jet, and (v) numerous enough to agree with the blazar flare duty cycle. 

Neglecting redshift factors in these estimates, we see that condition (i) implies that rcioud/cF^ ^ tf, 
so that the cloud radius rdoud <^ 3 x lO^^F^Qiday cm, where Fio = F/IO. Assuming that the plasmoid 
is highly radiative, equation (8) gives 4:Tt fbdl^^^^jUpT'^ cudoud ^ 10'**L48/fc, where dcioud = lO^'^rfie cm is 
the distance of the cloud from the supermassive black hole. Condition (ii) therefore implies that ricioud 
^ 2 X 10^L48(rio(ii6)~^ cm~^. In the fully radiative limit, we see from equation (20) that a plasmoid 
will be decelerated to nonrelativistic speeds if the number of swept-up particles roughly equals the 
number of jet particles. Consequently we find that the cloud mass Mdoud — -^plasmoid — Lftf/T, so that 
Mc\ond ^ W^'^L4sfbtday/Tio ergs. For a spherical cloud, we therefore see that ndoud ^ 10^L48/6iday/-Ri5Fio 
cm~"^, where we define R15 = rcioud/(10^^cm). 

For condition (iii), we find that the cloud takes ~ l-^dj^g /Mg' days to travel one degree if the clouds 
follow circular orbits, and could therefore traverse the jet opening angle on a time scale of several days for 
the given parameters. Condition (iv) is satisfied if i?i5 ^ diQ(9j/Q.l rad), where 9j is the jet opening angle. 
The duty cycle of fiares is not well known due to the sparse sampling of blazars, but might be in the range 
of 1-10%. Condition (v) therefore translates into a covering factor of comparable value as the duty cycle. 

The interpretation that blazar flares could result from relativistic plasma outflows interacting with 
clouds is therefore consistent with observations if marginally Thomson thick clouds of radii ^ 10^^ cm are 
located some hundreds to thousands of Schwarzschild radii from the central supermassive black holes. 



9. Summary 

A method is presented in this paper to treat the dynamics of a volume of relativistically moving 
magnetized plasma - a plasmoid ~ that sweeps up material from the interstellar medium and whose 
entrained nonthermal particles radiate their internal kinetic energy, thereby changing the relativistic inertia 
of the plasmoid. The formalism is cast in terms of particle distribution functions in the comoving plasmoid 
frame: these functions evolve when particles are swept-up, when energy is transferred between the hadronic 
and leptonic components in the plasmoid, and when energy is lost to radiation. Such an approach is 
open-ended in the sense that more complicated scenarios not dealt with here can be treated, such as the 
inclusion of other radiation processes (e.g., Compton scattering and bremsstrahlung, and secondary and 
photo-meson production), of diffusive escape of nonthermal particles from the plasmoid, and of pair and 
compactness effects and the microphysical details of particle acceleration within the plasmoid. 

We have studied a highly idealized system consisting of a relativistic blast wave or plasmoid with a 
constant entrained magnetic field. Only optically-thin synchrotron processes were treated here, and it was 
assumed that a fixed fraction of the swept-up nonthermal proton energy is instantaneously transformed 
into a power-law nonthermal electron momentum spectrum. This model was completely solved in the limit 
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of a rclativistic plasinoid whose Lorentz factor follows the power-law behavior given by equation (13). This 
solution can be used to benchmark numerical codes that include additional processes. Application of these 
results to studies of gamma-ray bursts and blazars was indicated, and will be developed more thoroughly in 
a numerical examination of the system. 

We thank M. Bottcher for criticisms which helped clarified the derivation of equation (6), and we 
acknowledge useful comments by M. G. Baring, A. K. Harding, and F. C. Jones. The work of CDD was 
supported by the Oflnce of Naval Research. The work of JC was performed while he held a National 
Research Council - Naval Research Laboratory Associateship. 

Table 1: List of Parameters, their Standard Values, and Values of Derived Quantities 



Parameter or 


Standard Value 


Derived Quantity 




r = ro(a:/a:o)"'' 


g = 3 


To = SOOraoo 


raoo = 1 


next = no(a;/a;())"'' 


r; = 


no 


1 cm^^ 


Eso = fio/lO"" ergs 


£50 = 1 


A = Aoix/xoY 


i = 2 


h 


1 


Ho 


10 Gauss 


s 


2 


z 


1 


gr ^j + 'i--v 


3 


9a = U + l- V)/2 


3/2 


Xo 


0.0018 pc 


Ao = 4:nxlft 


3.9 X 10^^ cm-2 


LO-^ 


156 sec 


fi-i 


0.58 sec 



Table 2: 


Negative of Slop 


es of Time Profiles: 


i in Fully Radiative Regi 


ime g : 


= 5r == 


J + 1- 


rj and 


Fully 


Non-Radiative Regi 


me g = 


ga = g, 


72 
























s = 1 








s = 1.5 






s = 2 






s = 3 




g 


1 


2 


3 




1 


2 


3 


1 


2 


3 


1 


2 


3 


Cl 


2/3 


3/5 


4/7 




2/3 


3/5 


4/7 


2/3 


3/5 


4/7 


2/3 


3/5 


4/7 


C2 


1/3 


2/5 


3/7 




1/3 


2/5 


3/7 


1/3 


2/5 


3/7 


1/3 


2/5 


3/7 


rad 
^uncooled 


1/3 


2/5 


3/7 




5/12 


1/2 


15/28 


1/2 


3/5 


9/14 


2/3 


4/5 


6/7 


rad 
•--cooled 


7/6 


6/5 


17/14 




5/4 


13/10 


37/28 


4/3 


7/5 


10/7 


3/2 


8/5 


23/14 


adi 

•-"uncooled 













1/12 


1/10 


3/28 


1/6 


1/5 


3/14 


1/3 


2/5 


3/7 


^cooled 


5/6 


4/5 


11/14 




11/12 


9/10 


25/28 


1 


1 


1 


7/6 


6/5 


17/14 


Ac 


5/6 


4/5 


11/14 




5/6 


4/5 


11/14 


5/6 


4/5 


11/14 


5/6 


4/5 


11/14 
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Explosion Frame 




A 

z 



V(x)^A(x)L(x) 



V(x+5x) = 
A(x+6x)L(x+dx) 



X 




x+5x 



%(><) 



A 
X 



r(x) 




r(x+8x) 



Fig. 1. — Cartoon illustrating the geometry of the plasmoid outflow and deceleration. The variation of the 
plasnioid bulk Lorentz factor T{x) with location x depends on the amount of swept-up external matter, 
which is characterized by density rioxt (x) , and the amount of internal energy which is radiated. 
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Fig. 2. — Solid and dotted curves show numerical solutions of equation (6) for the plasmoid Lorcntz factor 
r as a function of location x in the non-radiative and radiative regimes, respectively. Standard parameter 
values, shown in Table 1, are used. Solid and dotted lines give the power-law asymptotes in the respective 
regimes. 
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Fig. 3. — Dependence of F/Fo anid the ratio of distance traveled to distance xo given by equation (15) on 
the dimensionless quantity rjtobsj where 51 is given through equation (40) and iobs is the observer time. 
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Fig. 4. — Spectra (left) and time profiles (right) of a blast wave decelerating in a uniform external medium. 
Here it is assumed that the power- law behavior of the blast wave is described hy g ~ gr ~ 3, and that 
100% of the swept- up nonthermal proton energy is transformed into an equally energetic power-law electron 
distribution with a low-energy cutoff pmin = 0.2 and a high-energy cutoff pmax = 3 x 10^. All other parameters 
are standard (see Table 1). Observing times of the spectra are, from top to bottom, 1 sec, 10 sec, 100 sec, 
etc. The time profiles at 1 MeV, 1 keV, 1 eV, and 2.4 GHz are labeled by triangles, crosses, circles, and 
filled diamonds, respectively. 
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Fig. 5. — Same as Fig. 4, except that Hq = 1 Gauss. 
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Fig. 6. — Same as Fig. 4, except that Fq = 30, no = 0.01 cm '^, and fb = 0.01. 



